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We present further data in search of a parity-flavor-breaking phase in simulations of dynamical QCD with two 
flavors of light Wilson fermions in the strong coupling region. This is done on lattice sizes of 4 an up to 10* for 
a variety of values of (3 and k as well as the coefficient, h, of an explicit breaking term included in the action. We 
confirm the existence of a region in the (3 — n plane where such a phase exists at /3 = &/g^ smaller than 5.0 and 
above Kc- 



1. Introduction 

We extend our search for a parity-flavor- 
breaking phase in simulations of dynamical QCD 
with two flavors of light Wilson fermions to the 
strong coupling region of (3 less than 5.0. 

The signature used in this search is the same 
used in our previous work . We introduce an 
explicit breaking term into the action with a small 
coefiicient and then study the flnite size behavior 
of its expectation value as this coefficient tends 
to zero. 



2. Numerical Simulations 

We report here on simulations done with two 
flavors of Wilson fermions at (3 — 3.0, 3.5, 3.7 
and 4.0 on volumes of 4"^, and up to 10^ for a 
variety of values of k ranging from less than the 
appropriate Kc to values greater than k^. Some of 
these simulations (denoted by the letter C in the 
tables) are continuing and simulations a,t [3 — 4.5 
are in progress. 

We introduce into the QCD action a term of 
the form ihijj^^T^ip where ra is a 2 x 2 matrix 
representing the third element of the genera- 
tors of flavor SU(2) algebra. Upon integrating 
the fermionic variables this is reflected in the 



simulation by the product of two determinants: 
DetM(h) * DetM{-h) where M (h) is given by a 
simple modification of the Wilson Matrix M^, as: 

M{h) = + ih-f5 

As pointed out by Aoki ||], we also have here 

75M(-/i)75 AP 



and: 



DetM{-h) = DetM\h) , 
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Simulations were done for the parameter h 
taking values ranging from 0.005 to 0.1. For the 
volume dependence we concentrate on the smaller 
values of h and in particular h = 0.005 for all 
three volumes considered and mostly for values 
of K greater than Kc. 

The order parameter we compute is the expec- 
tation value of the operator i^^^r^ip. With our 
notation this is given as 

PF = -/mTr(75M-i(/i)) 

3. Results 

For the values of /3 considered, simulations were 
performed, as mentioned above, at various values 
of K both below and above Kc- We shall present 
the data and results for each value of f3 considered 
separately. The results are shown in Tables 1-4. 
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Table 1 

Parameters and measured order parameter PF^ 
for the case oi (3 = 3.0 on lattices of volume L**, 
for L = 4,6, and 10. 



K 


h 






PFu) 


0.1900 


0.005 




0.159(3 ) 


0.159(1) 


0.1900 


0.050 




1.503(24) 


1.503(8) 


0.1900 


0.100 




2.670(37) 




0.220 


0.005 




0.54 (4) 


0.526 (14) 


0.220 


0.050 




2.681(73) 


2.690 (26) 


0.220 


0.100 




3.654(64) 




0.250 


0.005 


2.13 (64) 


2.60 (30) 


2.99 (13) C 


0.250 


0.050 




3.67 (11) 




0.250 


0.100 




4.232(76) 





Table 2 

Parameters and measured order parameter PF^ 
for the case of /? = 3.5 on lattices of volume 
for L = 4, 6, 10. 





h 


PFi 


PFe 


PFw 


0.1900 


0.005 




0.1648(33) 


0.1639(11) 


0.1900 


0.050 




1.531 (26) 


1.533(9) 


0.1900 


0.100 




2.682 (38) 




0.240 


0.005 


1.42 (60) 


2.16(28)0 


2.43(9)C 


0.240 


0.050 




3.171(94) 




0.240 


0.100 




3.848(79) 




0.250 


0.005 


0.68(42) 


0.354(5)C 





It is clear that if one concentrates on the behav- 
ior with volume of the simulations aX h = 0.005, 
one may quickly note that a signature for the 
parity-flavor-breaking phase is indicated for the 
following values of /3 and n: (3.0, 0.25), (3.5, 
0.24), (3.7, 0.24) and (4.0, 0.24). See Fig. |l| for 
the latter two points. 

This is also confirmed by following the behavior 
of the order parameter as a function of k at 
constant volume and constant (3. it is clear that 
all values above the order parameter goes through 
a relative maximum. 

The behavior of the order parameter at all 
other values of (3 and k, computed in Tables 1-4, 
are consistent with an absence for a signature of 



Table 3 

Parameters and measured order parameter PF^ 
for the case of /3 = 3.7 on Lattices of Volume L^, 
with L ^ 4,6,8, 10. 



K, h PFi PFe PFio 

0.2300 0.005 1.14 C(45) 1.78(26)0 1.8(2)C 
0.2400 0.005 0.673(380) 1.128(239) 2.43(10)C 

PFa 

0.2500 0.005 0.337(193) 0.360(62) 0.355(63T 



Table 4 

Parameters and measured order parameter PF^ 
for the case oi (3 — 4.0 on Lattices of Volume L"*, 
with L = 4,6, 10. 



K 


h 


PFi 


PFe 


PFio 


0.1900 


0.005 




0.1715 (38) 


0.1719(13) 


0.1900 


0.050 




1.574 (31) 


1.569 (10) 


0.1900 


0.100 




2.691 (38) 




0.2200 


0.005 


0.58(21) 


1.20 (18) 


1.25(07)C 


0.2200 


0.050 




2.605 (76) 




0.2200 


0.100 




3.358 (64) 




0.2300 


0.005 


0.18(03) 


0.33C(11) 




0.2400 


0.005 


0.28(13) 


1.30(23)C 


2.38(10)C 


0.2400 


0.050 




2.459 (86) 




0.2400 


0.100 




3.517 (68) 





a parity-flavor-breaking phase. This behavior is 
in fact very similar to the behavior for values of 
/3 > 5.0 which were discussed in ref. jl|. 

It is important to note here that the phase 
is not indicated for some points at the same /3 
but larger values of k. In particular the phase 
disappears at (3.5, 0.25), at (3.7, 0.25) . 

It is clear from the discussion above that QCD 
with two flavors of Wilson fermions does exhibit 
a parity-flavor-breaking phase at values of (3 up 
to 4.0 , as postulated by Aoki and collaborators. 
It has also been demonstrated in [|] that this 
phase does not seem to extend beyond (3 = 5.0 . 
This information is exhibited in Fig. ^. It is then 
possible that this phase pinches out, in a manner 
similar to that in the NJL model, at < 5.0. 
In this case, this phase would not be relevant for 
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^= 3.7 K= 0.24 
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Figure 1. Histogram signature of a parity-flavor- 
breaking phase for /? — 3.7 and 4.0; k — 0.24. 



the discussion of the approach to the chiral Umit 
in QCD and the ensuing Goldstone nature of the 
pions for /3 > 5.0. In fact, aU indications are 
such that, as shown formally sometime ago, this 
is simply related to the approach to zero lattice 
spacing and infinite volume. 
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Figure 2. Parity- flavor phase diagram. Continu- 
ous line represents approximate positions for chi- 
ral and finite temperature transitions at TVj = 4 
for comparison. 
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